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The gravitational effect on tfie chiral phase structure at finite temperature and density is in- 
vestigated in the framework of Nambu-Jona-Lasinio model in D-dimensional ultrastatic Einstein 
universe. In mean field approximation, the thermodynamic potential and the gap equation deter- 
mining the curvature, temperature and density dependence of the chiral condensate are analytically 
derived. In the sense of chiral symmetry restoration and the order of the phase transition, the scalar 
curvature of the space-time plays similar role as the temperature. 
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I. INTRODUCTION 

The phase structure of quantum field theory, espe- 
cially quantum chromodynamics(QCD), was widely in- 
vestigated in fiat space-time in the last two decades. It 
is quite interesting to make studies of how the phase 
structure changes in the circumstance of very compact 
stars and early universe where not only the temperature 
and density effect but also the gravitational and gauge 
field effects can not be neglected. The QCD phase tran- 
sitions that probably happen in these systems include 
mainly the restoration of spontaneously broken chiral 
symmetry and the deconfinement from hadron gas to 
quark matter and vice versa. They are closely related 
to the hadrosynthesis at the early universe. We focus 
on the chiral phase transition in this paper. 

Four fermion models (and gauged four fermion mod- 
els) are often applied to describe the dynamical sym- 
metry breaking in QCD, electroweak theory and grand 
unified theory. One of the most frequently used 
four fermion models is the Nambu-Jona-Lasinio (NJL) 
modelpy]. The chiral symmetry of this model is sponta- 
neously broken according to the emergence of a vacuum 
condensate of the composite field of fermions, {"tpi/j) ^ 0, 
and the fermion mass is proportional to the conden- 
sate. The NJL model is particularly convenient for the 
investigation of chiral symmetry when some external 
conditions, like temperature, chemical potential, non- 
trivial topology, external gauge field and gravity field 
are taken into account. Such conditions need to be con- 
sidered in the study of high energy heavy ion collisions 
where the temperature and baryon density are high 
enough, compact stars like neutron stars and quark 
stars where the baryon density may be ten times the 
normal nuclear matter density and the strong gravity 
and external magnetic field can not be neglected, and 
early universe where the gravity is certainly important 
and the temperature is extremely high. Because of the 
non-perturbative phenomena around the phase transi- 
tion and the fermion sign problem in lattice QCD at fi- 
nite chemical potential, the NJL model has been widely 
considered as a low energy effective theory of QCD in 
the study of chiral phase transition in flat space-time [^. 

The chiral symmetry breaking in four fermion models 
in curved space-time at zero temperature and chemical 
potential was investigated in recent years [1, 0, 0, 0, 
0, H, H, • It was found that when the curvature of 
space-time is sufficiently large, there will be no more 
spontaneous chiral symmetry breaking in the vacuum. 



Therefore, in the sense of chiral symmetry restoration, 
the curvature R behaves like temperature T. This phe- 
nomenon can be clearly seen in the discussion with both 
curvature and temperature effects [I. [ill, [l^. With the 
method of weak curvature expansion for fermion prop- 
agator, the chiral phase structure of Gross-Neveu|13l| 
model in 2-t-l dimensional space-time is studied at fi- 
nite chemical potential /i but zero temperature . the 
phase diagram in the R — fj, plane is similar to the dia- 
gram in the T — ji plane in flat space-time. In the four 
dimensional curved NJL model at finite temperature 
and chemical potential, the weak curvature approxima- 
tion leads to a first-order chiral phase transition[15]. 

A natural question we ask is if we can use the method 
of weak curvature expansion to treat the chiral prop- 
erties around the critical curvature Rc, when Rc is not 
small. In this paper, we study how the temperature, 
chemical potential and external gravity influence the 
chiral symmetry of the NJL model in different dimen- 
sions. We take the Einstein universe as the background 
space-time, in which we can analytically calculate the 
thermodynamic potential without using the weak cur- 
vature approximation in mean field approximation. In 
flat space-time, the chemical potential effect is very dif- 
ferent from the temperature effect on the chiral phase 
transition. For instance, the transition is of first order 
at high chemical potential, while it is of second order at 
high temperature. We want to know how the chemical 
potential affects the chiral phase transition in highly 
curved space-time. 

The paper is organized as follows. In Section [III 
we review the NJL model in curved space-time at fi- 
nite temperature and chemical potential and obtain the 
thermodynamic potential in mean field approximation. 
In Section imi we analytically derive the thermodynamic 
potential and gap equation for the chiral condensate in 
the Einstein universe. The phase diagrams in differ- 
ent planes are shown and discussed in Section HVl We 
summarize in Section fVl 



II. THE CURVED NJL MODEL IN MEAN 
FIELD APPROXIMATION 



In _D-dimesional curved space-time with signature 
(+,-,-,-,•■ •), the action S of the SU{N) NJL 



2 



model takes the form of 



G / - 



(1) 

where g is the determinant of the metric tensor g^^, (x) , 
ipix) represents a fermion field in a A'^-dimensional in- 
ternal space, "i/j — ■0^7" is the corresponding Pauli con- 
jugate spinor, the index with a hat denotes a standard 
Dirac matrix defined in Minkowskian space-time, G is 
the four fermion coupling constant, j^ix) is the Dirac 
matrix in curved space-time satisfying {7^^,71/} = gfj.^, 
Vjy is the covariant derivative with the chemical poten- 
tial parameter /i, Vjy = — i^Soi, = + — iuSqi, 
with riy{x) being the spinor connection [TgI, [20I [21I. [23| . 
and 7'^ is defined as the same as in flat space-time, 
it is independent of the space-time coordinates. The 
model has a global SU{N) (g) U{1)l U{1)r symmetry 
at /X = 0. 

It is convenient to introduce the tetrads 
formulism [H, [2l|, to connect the defini- 

tions in curved and flat space-times. Let e°(a:) be the 

tetrads defined by = Vab^li^t' where 77.^ is the 
Minkowskian metric, we can rewrite 7'^ = 7^6^^, and 

with = d^el + r-^e? 

0, where e'^ is 



fib 

V^7' 



under the requirement V^7'' 
the inverse of e° . 

The essential quantity characterizing a system in a 
grand canonical ensemble can be taken to be the par- 
tition function Z. It can be expressed in terms of the 
action of the system as 



(2) 



After bosonization, it becomes the integration over the 
boson fields a = —2G'ipip / N and tt = —2G'ipij5ip / N , 



Z = J [da] [dir] , 



with the effective action 



N 



(3) 



-i In Det 



4G 

n^V^ - (cr + ^7^71-) , (4) 



where the determinant is taken in the internal space, 
the Dirac space and the space-time manifold. 

In flat space-time, in order to incorporate finite tem- 
perature effects in quantum field theory, we can per- 
form a Wick rotation t —> —it and J dt —>■ —i dr 
with (3 = 1/T, and the integration over fields is con- 
strained so that "0(0) = — and ^(0) — — 
are satisfied for fermions. This is the so-called imag- 
inary time formalism of finite temperature field theory. 
However, in curved space-time with non-static metric, 
the presence of time-space cross terms makes it diffi- 
cult to have a well-defined Wick rotation. Only for the 
ultrastatic metric (for the static metric, we can con- 
formally transform it to the case of ultrastatic metric, 
see refs.fl^ 113), we can have a well-defined temper- 
ature by a Wick rotation. The parameter has a 



clearly physical meaning of chemical potential, since 
the number density is still conserved in curved space- 
times, d/dt J d^~^a;-\/|g|^7°?/; — 0. We will discuss 
only the ultrastatic universe M ® E below with E being 
the {D — 1) dimensional space manifold. 

With the known partition function, we obtain the 
thermodynamic potential 51 as a function of tempera- 
ture and chemical potential. 



n{T, n) 



T 

-— lnZ(r,^), 



(5) 



where V ~ J d^^^x^/\g\ is the volume of the space 
manifold E. 

In mean field approximation which is just the 0(1) 
order in the large TV expansion of fermion propagator, 
the fields a and tt are replaced by their thermal aver- 
ages (cr) and (tt). To simplify the notations, we express 
in the following the averages (a) and (tt) by a and tt, re- 
spectively, without making confusion in the mean field 
approximation. With the short notations, the thermo- 
dynamic function can be simplified as 



AG 



T 

TT^) - — In Det 



«7^Vp - {<J + i7^7r) 



(6) 

If we do not consider the chiral anomaly induced by 
gravity, the thermodynamic potential is invariant under 
chiral transformation, and we can set tt = without 
loss of generality. In this case, we have 



Na^ T 



AG 



V 



■In Det 



Z7^V, 



(7) 



The thermal expectation value a is the order param- 
eter of the chiral phase transition. The thermodynamic 
potential ([7]) is a function of T and /i with a initially 
an undetermined parameter. In the spirit of thermo- 
dynamics, the physical system is described only by T 
and /I. The order parameter as a function of T and /i 
is determined by the minimum of the thermodynamic 
potential^, i26i] . 



do. 



0, 



9^ 



> 0. 



(8) 



Since we did not consider the initial fermion mass in the 
action ([1]) , the effective fermion mass produced through 
the spontaneous chiral symmetry breaking is just the 
chiral condensate, to = cr. 

The measurable bulk quantities like pressure p, en- 
tropy density s, fermion number density n, and energy 
density e are related to il by 



p = —Q, s = 

e = —p + Ts 



_dn 



(9) 



Now we come back to the calculation of the mean 
field thermodynamic potential ([7]). In ultrastatic space- 
time with line element ds^ — dt^ — gij{x)dx^dx^ , it is 
easy to show that the zero component Fq of the spinor 
connection is zero, and we have g^'^V fiVi, = d"^ jdt^ — 
with V'^ = g^^^ iS j being the spinor Laplacian de- 
fined in the space manifold E. Taking into account the 
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relation -if^Y^t^y = s'^'^V^V^ + \R- - 2i/iVo, 
we have 



21nDet 
InDet 
InDet 

InDet 



j7"V^ - a 
lY^^ - a 



InDet 



d_ 

dt 



V^ + ^ + a' 
4 



(10) 



In order to complete the calculation, we should solve 
the eigenequation of the operator in the last square 
bracket, 



d_ 

di 



111 
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(11) 



where k stands for a set of complete quantum numbers 
describing the eigenvalue A and eigenstate 5'. By sep- 
arating the t and x dependence of the x dependent 
part 0/ (x) satisfies the corresponding eigenequation 



R 



t>i = Ui (pi 



(12) 



in the space manifold and the normalization condition 

d''-^xy^\M^)cl>i'{S)=5n'. (13) 

If the solutions ujf and (l)i{x) of the eigenequation 
(|12p are known, we may take 



Afe ujf + (J^ - (po - lif, 
^k{t.x)^e'^^«'4>m. 



(14) 



where po = «(2n + l)7r//3 is the fermion frequency with 
n = 0,1,2,..., since the eigenfunction must be an 
antiperiodic function of the period 13. 

After the frequency summation over n, we finally 
obtain the thermodynamic potential in mean field 
approximation [H, [2^ H^l 



AG 
+Tln(l 



+Tln 1 



(15) 



where Ei = ^ u]f + is the quasiparticle energy and 
di is the degeneracy of the Z-th eigenvalue of equation 

(ini. 



III. THE GAP EQUATION IN THE EINSTEIN 
UNIVERSE 

The question left is to solve the eigenequation ([T^ . 
To this end, some approximate methods like weak 



field approximation!^ [^, [TJ, [T^ and high temperature 
expansionlil, [ll] are widely used. However, around 
the phase transition, the interaction among the con- 
stituents of the system should be very strong and those 
methods based on perturbative expansion are in prin- 
ciple not suitable for the study of the phase structure. 
In fact, the eigenequation can be exactly solved 
in some special universes ^5, 28, 29], for example, the 
static Einstein universe. The Z^-dimensional Einstein 
universe is represented by the line element 



dt^ -a^{d0^ +sm^ ednD-2 



(16) 



defined with the topology R ® ^, where a is re- 
lated to the curvature, R = {D - 1){D - 2)a~^. With 
the eigenvalues ^28, 29] of the spinor Laplacian on the 
sphere §^~^, the degeneracy di, the eigenvalue and 
the volume V can be expressed in terms of I and R as, 

2[(^+^)/^ir(/ + iJ-i) 
^ nT{D - 1) ' ' ^ ' ~ ' 

= + ^ = 0,1,2,.., 

' 4(1) - \){D - 2) 4' ' ' ' ' 

, , 2^^/2a^-i 



(17) 



where \x\ is the floor function of a;, and the summation 
over the eigenvalues in the mean field thermodynamic 
potential becomes now explicit. 



r2 = — — d]e 

4G V ^ ^ 



+Tln 1 



-Tin 1 



(18) 



Because of the contact interaction among the particles, 
the NJL model is non-renormalizable when D > 4, and 
in general it is necessary to introduce a regulator that 
serves as a length scale in the problem. That is the rea- 
son why we have introduced a soft cutoff factor e~^>'l ^ 
in the summation over I. Otherwise, the first term in 
the square bracket which is the contribution from the 
vacuum will be divergent. 

There are two parameters in the model, the coupling 
constant G and the soft cutoff A. In the case of flat 
space-time, they can be fixed by fitting some observ- 
able quantities in the vacuum, for instance, the pion 
mass and the pion decay constant. Since we are not 
familiar with the particle properties in curved spaces, 
we are in principle not able to determine the two pa- 
rameters in the Einstein universe. However, what we 
want to do in this paper is a general study of the grav- 
itational effect on the chiral phase structure, without 
considering a specific system. We focus only on the 
qualitative dependence of the phase structure on the 
gravity, and do not care very much the precise criti- 
cal values of the curvature, temperature, and chemical 
potential. To this end we scale the thermodynamic po- 
tential (UHl) by the soft cutoff A, 



4 



n _ N{a/Kf N 



1=0 



El T, , 



Ei/A-,i/A 

■ e 



y / E,/A + n/A 

- In I 1 + e ^/'^ 



(19) 



Without causing confusion, we represent the scaled dimensionless quantities il/A^, cr/A, A^^^G, A^^^V, 
R/A'^, T/A, n/A, oJi/A still by the corresponding quantities il, a, G, V, R, T, ji, oji, and rewrite p9|) in 
the dimensionless form, 



4G V 



Ei+Tlnil 



1=0 



(l + e — \ +rinn + e —\ 



(20) 



r 



Now only one dimensionless parameter, the coupling 
constant G appears in the model. 

Calculating the first order derivative of the thermo- 
dynamic potential ((20|) with respect to the condensate 
(T, we obtain the gap equation which determines the 
T— , \i— and i?— dependence of the condensate, 



cr[l - 2GI{T,ii,R,a)\ = 0, 



(21) 



with the function / defined as 



V El 

1=0 



{l-f{Ei+fi)-fiEi-fi)), (22) 



where f{x) — l/(e^/^ + 1) is the Fermi-Dirac distri- 
bution function. The trivial solution cr = of the gap 
equation describes the symmetry restoration phase, and 
the other solution a ^ corresponds to the symmetry 



breaking phase. It is necessary to note that the solution 
of the gap equation is not guaranteed to be the physical 
condensate, we should check for each solution if it is the 
position of the minimum thermodynamic potential. 

What happens when the curvature tends to zero? 
While the line element (fT6)) of the Einstein space-time 
can not return to the line element of Minkowskian 
space-time when — > or a — > oo, because the two 
space-times have different topologies, the thermody- 
namic potential in the limit R ^ becomes the one 
in the flat space-time. In the limit of a oo, we have 



di 



a ^l, duji = uji+i - uji 

2l(D+l)/2]^D-2^D-2 



(23) 



and the thermodynamic potential is reduced to 



J 



n 



Na^ iV2[(^-i)/2lr(£)/2) 



4G 



7r^/2r(D- 1) 



dpp 



£;„ + Tln 1 



Tin 1-Fe" 



(24) 



where we have replaced the integrated variable wi by p 
which could be considered as the fermion momentum. 
This is exactly the familiar formula for the thermody- 
namic potential of the NJL model in Z3-dimensional 
Minkowskian space-time. 

Before discussing the chiral properties at finite tem- 
perature and density, we first consider the pure gravi- 
tational effect. In the curved vacuum with T = /i = 0, 
the thermodynamic potential is reduced to 



4G V 



Y.diE. 



le 



(25) 



1=0 



and correspondingly, the gap equation becomes 



= fj 



2G^ di 



(26) 



To have chiral symmetry breaking a ^ in the vacuum, 
the coupling constant G must exceed the critical value 



Gc which is determined by 



_ 2Gc di _ 



(27) 



The critical coupling constant Gc as a function of the 
curvature R is shown in FiglT] for the dimension D = 

4, 5 and 10. As i? increases, Gc growths monotonously, 
reflecting the fact that the broken chiral symmetry can 
be restored by finite curvature effect. To keep chiral 
symmetry breaking in curved vacuum with finite R, it 
needs stronger interaction than that in flat vacuum. 
However, this curvature effect is gradually washed out 
when the dimension of the space-time increases. From 
FiglU the effect is almost fully cancelled in the case 
of I? = 10. To have a remarkable curvature effect, we 
consider in the following only the cases with D = 4 and 

5. The initial value of Gc at i? = can be analytically 
derived by taking the limit a ^ oo in Ea. (P7|) . 



Gc(i? = 0) = 



(i:>-2)7r^/22-[(-0+i)/2l 

r(^/2) ■ 



(28) 
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FIG. 1: The critical coupling constant Gc as a function of R 
for three values of D. At each value of D the chiral breaking 
phase corresponds to the region above the line. 
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FIG. 2: The shifted thermodynamic potential as a function 
of chiral condensate at T = = for 7? = 0, 2, 10 at D = 4 
and 7? = 0, 2, 12 at D = 5 



IV. PHASE DIAGRAMS 



To guarantee chiral symmetry breaking in the vac- 
uum, we take the coupling constant G = 10 > Gc(0) 
for D — A and 5 in the following numerical calculations. 

The curvature effect on the chiral symmetry restora- 
tion can be clearly seen in the i?— dependence of the 
thermodynamic potential. In Fig[2]we plot the shifted 
thermodynamic potential in the curved vacuum, U,{a) — 
il(0), as a function of a for different values of R. At 
small values of R, the physical condensate correspond- 
ing to the minimum thermodynamic potential is finite 
but not zero, which means chiral symmetry breaking, 
but at large values of i?, the minimum of the thermo- 
dynamic potential is located at ct = 0, which stands 
for chiral symmetry restoration. It is easy to see from 
Figs[T] and [2] that the critical curvature Rc for chiral 
restoration increases with increasing dimension. 

We now calculate the chiral condensate at finite tem- 
perature, chemical potential and curvature. At T = 0, 



the gap equation (|2T|) is reduced to 



1 



'V ^ El 

1=0 



= 0, 



(29) 



where 9{x) is the step function. Since Ei grows 
monotonously with increasing I, the nonzero chiral con- 
densate is just a constant ctq in the region 



M < fJ-o = Eq 
with do determined by 



(=0 



0. 



(30) 



(31) 



However, as we mentioned above, the solution /ig of the 
gap equation is probably not the critical chemical po- 
tential /ic for chiral phase transition. We should check 
the minimum of the thermodynamic potential 



2 oo 

n{p, i?, a) = ^ [/i + {El - e{Ei - /i)] 

(32) 



1=0 



and determine by the condition 

r!(/i„i?,(To) = r!(/x„i?,0). 



(33) 



At i? = 0, we have = 2.2 < /iq 2.9 for L» = 4 
and /ic — 3.5 < fj,o = 4.7 for D — 5. At fixed R, the 
physical condensate a is simply a step function of /i. 



cr(/i) = (ToO{^j.c " ^J■)■ 



(34) 



We show the chemical potential dependence of a in 
FigElat r = and for i? = and 3. The critical chemi- 
cal potential drops down with increasing curvature and 
finally reaches y^ujQ{Rc) + (Tq{Rc) at the critical curva- 
ture Rc- It is easy to see that Rc satisfies the condition 



V ^ LUl 

1=0 



0. 



(35) 



At D = 4 and 5, the critical curvature Rc are, respec- 
tively, 4.7 and 8.5. Above Rc there exists no more chiral 
breaking phase and the condensate keeps zero at any 
chemical potential. 

The chiral condensate at finite temperature is shown 
in Fig|4] at fixed R. Similar to the case in flat space- 
time, the temperature effect leads to chiral symmetry 
restoration. When /i = 0, the chiral condensate de- 
creases with increasing temperature and finally reaches 
zero at the critical value Tc determined by 



.G^di 
uji 



{l-2f{uji)) = 0. 



(36) 



Above Tc the chiral breaking phase disappears at any 
chemical potential. 

The phase diagram in the T — fi plane at fixed R is 
shown in FiglH It is very similar to the familiar phase 
structure in fiat space-time. Each phase transition line 
separates the chiral breaking phase below the line and 
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FIG. 3: The chiral condensate a as a function of chemical 
potential /i at T = and for i? = 0, 3 and D = 4, 5. 
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FIG. 4: The chiral condensate a as a function of chemical 
potential ^ at fixed curvature R = 1 and for T = 0, 1, 1.3 at 



the chiral restoration phase above the line. The chiral 
breaking region is gradually suppressed with increasing 
curvature and finally disappears at the critical value 
Re- The dashed and solid lines represent, respectively, 
the second and first order chiral phase transitions, and 
the dots indicate the tricritical points which link the 
corresponding second and first order transitions. 

The phase diagrams in i? — /i and R — T planes are 
shown, respectively, in Figs [5] and [71 Again the dashed 
and solid lines mean the second and first order phase 
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FIG. 5: The phase diagram inT — ^ plane for R = Q,2, 4.5 
at D = 4 and R = 0, 4.5, 8 at D = 5. The dashed and 
solid lines represent, respectively, the second and first or- 
der phase transitions, and the dots indicate the tricritical 
points which connect the corresponding second and first or- 
der transitions. 
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FIG. 6: The phase diagram in i? — ^ plane for T = 0, 0.5, 1 
at 7? = 4 and T = 0, 1, 2 at D = 5. The dashed and sohd 
lines represent, respectively, the second and first order phase 
transitions, and the dots indicate the tricritical points. 



transitions, and the dots indicate the tricritical points. 
Just as what we expected, at fixed curvature, the phase 



7 




D=5 




0.5 1 1.5 2 2.5 3 



FIG. 7: The phase diagram in R — T plane for /i = 0, 1.7, 2 
at D = 4 and /x = 0, 2, 3 at D = 5. The dashed and solid 
lines represent, respectively, the second and first order phase 
transitions, and the dots indicate the tricritical points. 



transition becomes more and more easy when the tem- 
perature or chemical potential increases. The straight 
dashed lines in the R ~ ji plane at T = reflect the 
fact that the Einstein universe is compact and in turn 
the lowest eigenvalue luq of equation (|12p is not zero. 
The toothed lines in FigHat T = come from the step 
function in equation (|29p which makes the low limit of 
the summation jump up when /i increases. From the 
phase diagrams in FigslSl [6] and [71 the second order 
phase transition happens at high T or high R or low ^, 
while the first order phase transition occurs at low T or 
low R or high fi. While finite values of R, T and /i can 
all lead to a chiral phase transition, the order behav- 
ior of the phase transition induced by curvature effect 
is similar to that by temperature effect, but different 
from that by chemical potential effect. Therefore, in 
the sense of chiral restoration, the gravitational effect 
is more like the temperature effect. 

We now discuss the critical properties for the sec- 
ond order chiral phase transition. Expanding the gap 
equation for a in the chiral breaking phase 



1 - 2GI{T,fi,R,a) = 
around the critical value Rc at fixed T and fi, 

R — Rr 



Rr 



b<T^ + CO--* = 



(37) 



(38) 



with the coefficients a, b and c defined as 



dl 

a{T,n,Rc) = 
b{T,fj,,Rc) = ^ 



a=0 



c(T, Rc) 



1 d^I 



2 9(f72)2 



(T = 



(39) 



at the critical point, where we have considered the fact 
that / is a function of ct^ and neglected the higher or- 
ders of cr^. 

The critical condensate CTc is then determined in the 
limit R Rr, 



= 0. 



(40) 



For a second order phase transition, there should be 
only one solution — 0, the coefficient c must vanish. 
Therefore, around the critical point, a behaves as 



a{T,ii,R) = ^- 



R — Rr. 



Rr 



1/2 



(41) 



which means the critical exponent Pn — 1/2. 

For a first order phase transition, however, there 
should be two values of at the critical point, one 
is CTc = 0, and the other is CTc 0. From the 
nonzero condensate is (Tc = ■\/— 6/c. Approaching to 
the tricritical point which is the end of the first order 
phase transition and the beginning point of the second 
order phase transition, the nonzero CTc — > 0, we have 
therefore 



6(T, /i, Rc) = 



(42) 



at the tricritical point. In this case, a around the tri- 
critical point behaves as 



<j(T,fi,R) 



R — Rr 



Rc 



1/4 



(43) 



which means = 1/4 at the tricritical point. Ob- 
viously, the discussion above for curvature is also valid 
for temperature, we have /3r = 1/2 for the second order 
phase transition and /3t = 1/4 at the tricritical point. 

The location (T^, /ic) of the tricritical point at fixed R 
is determined through its definition (|42p together with 
the gap equation. 



l^2GI{Tc,^ic,R,0)^0, 
b{Tc,^ic,R) = Q■ 



(44) 



FiglH shows the tricritical Tc and ^. c as functions of R . 
Tc ends at zero, while ^c ends at ujq(Rc) + <Jq{Rc). 



V. SUMMARY 

We have investigated the gravitational effects on the 
chiral properties at finite temperature and density in 
the framework of the NJL model in the Einstein uni- 
verse. In the mean field approximation, the thermody- 
namic potential and the gap equation to determine the 
chiral condensate are analytically derived without any 
further approximation, which allows us to reasonably 
study the chiral symmetry around the phase transition 
where any perturbative expansion is in principle not 
valid. 

If we take the curvature R of the space-time which 
describes the gravitational effect in the Einstein uni- 
verse as an external parameter, the role it plays in the 
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R 

FIG. 8: The tricritical temperature and chemical potential 
as functions of the curvature for D = 4 and 5. 



chiral phase transition is very much hke the effect of 
temperature. Both finite R and T result in chiral phase 
transitions and the behavior of the order of the transi- 



tion is similar. However, this gravitational effect on the 
phase transition will be washed out when the dimension 
of the universe is high enough. 

The chiral phase diagram in the T — /i plane is similar 
to the one in flat space-time where the thermodynamics 
can be obtained from that in the Einstein universe by 
taking the limit i? — > 0. However, the compact property 
of the Einstein universe makes the minimum critical 
chemical potential fic sd, T — he finite but not zero, 
which leads to a flat roof structure in the phase diagram 
in the R ~ fi plane at T = 0. From the definition of 
the tricritical point, we determined its location in the 
T — /i, R — fj, and R — T planes. 

It should be noted here that in the problem of the 
early universe the effects of external gauge field on the 
chiral phase transition are important. These effects 
have been studied in [13, [H, [3] , and our formulation 
permits the inclusion of these effects. 

The change in the chiral properties induced by finite 
curvature effects may be useful for the investigation of 
compact stars where the gravitational effects can not be 
neglected. A natural extension of the study is to discuss 
the di-fermion condensate in curved space-time, which 
is significant for the study of color superconductivity, 
since it may exist in the core of compact stars. 
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